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OUTLINE 

• National Institute of Aerospace (NIA) overview 

• NASA Langley (LaRC) overview 

• Examples of composites simulation: 

> Thermo-mechanical material model 

> Damage analyses of composites 

> Progressive damage material model 

> Virtual crack closure technique (VCCT) 

> Decohesion element 

> Flight 587 structures investigation 

> Rotorhub flexbeam analysis 

> Mixed-mode delamination failure criterion 

> Delamination in z-pin reinforced laminates 

• Concluding remarks 
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NATIONAL INSTITUTE OF AEROSPACE 


• An Independent Non-profit Research and Graduate Education Institute formed in 
2002 by a Consortium of Six Universities and the AIAA Foundation 

• Conceived by NASA Langley Research Center and established to serve as LaRC’s 
Collaborative Partner 

• Conducts Collaborative Research in Engineering and Science relevant to Aerospace 


Georgia Tech 
Virginia Tech 
University of Virginia 
University of Maryland 
North Carolina A&T State University 
Old Dominion University 
College of William & Mary 
Hampton University 
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COMPOSITE RESEARCH ACTIVITIES AT 
NASA LANGLEY 


• Computational materials 

• Crash worthiness of composite structures 

• Structural tailoring with composites 

• Manufacturing and fabrication technology 

• Health monitoring using embedded sensors 

• Residual strength and damage propagation 

• Influence of generalized imperfections on composite shell response 

• Delamination and crack growth 

• Thermo-mechanical material response 

• Multidisciplinary design environments 

• Uncertainty quantification for composite designs 

• Impact response and strain rate sensitivity 
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SIMULATION ISSUES FOR COMPOSITES 


• Variabilities in composite design (fiber placement, fiber angle, thickness, 
volume fraction, failure modes, etc) 

• Visualization of composite simulation results 

• Micro-mechanics through macro-mechanics - Problem of scale (global local 
analysis requirements) 

• Computational models for new and evolving materials 

• Computational models for incorporating mechanical-based failure models 

• Failure initiation and damage propagation of different composite architectures 
(sandwich construction, integrally-stiffened sections, etc) 

• Corroborating experimental program for validation of analysis 

At NASA Langley, these issues are addressed by researchers using: 

1 . User-defined material models 

2. User-defined element routines 

3. User-developed pre and post-processing software 
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NASA LANGLEY RESEARCH CENTER 


Research Technology Directorate (RTD) consists of 21 branches: 
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THERMO-MECHANICAL MATERIAL 

MODEL 


TRAVIS TURNER 


STRUCTURAL ACOUSTICS BRANCH 
NASA LANGLEY RESEARCH CENTER 

HAMPTON, VA 
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NUMERICAL MODEL 



Shell Element Mesh 



Pre-preg only 
• Plies 0.0036” thick 


Pre-preg & actuators 

• Plies 0.0032” thick 

• Actuators 0.006” thick 
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• Developed thermo-mechanical FE model based upon LaRC-developed 
constitutive model implemented in MSC.Nastran and ABAQUS 

• Shell element mesh separates glass-epoxy-only and SMAHC element types 


• Nonlinear static solution performed with imposed temperature load specified by 
experimental measurements at critical temperatures 
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STRUCTURAL DEFLECTION CONTROL 

Actuators embedded within layers of a laminated composite structure 


Machined shape 



Completed Flow Effector 




• Discrete actuator “inclusion” 

• Non-uniform consolidated pre-preg ply thickness 

• Non-uniform temperature distribution in service 
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Flow Effector Deflection Control 


Bench Top Test System 



Tip Deflection Comparison 
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Excellent numerical-experimental agreement 
Numerical design tool validated 
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DAMAGE ANALYSES OF COMPOSITES 



Through-the-thickness crack 
•fracture mechanics and modifications 
•strain softening 


Ply Damage 

•continuum damage modeling (CDM) 
•strength-based methods (criteria) 
•micromechanics approach 


Delamination/Debonding 

•fracture mechanics approaches (VCCT) 
•decohesion elements 
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PROGRESSIVE FAILURE MATERIAL 

MODEL 


NORMAN F. KNIGHT 

GENERAL DYNAMICS 
ADVANCED INFORMATION SYSTEMS 

CHANTILLY, VA 

Resident at Durability, Damage Tolerance and 
Reliability Branch, NASA Langley Research Center 
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PFA/PDA SIMULATIONS 


• Develop user-defined material subroutine UMAT for PFA and user-defined 
element subroutine UEL for PDA using ABAQ US/Standard 

• UMAT features linear elastic, bimodulus, orthotropic material model for a 
composite laminate 

• Failure initiation based on material allowable values using: 

- Maximum stress criteria 

- Maximum strain criteria 

- Hashin criteria 

- Tsai-Wu polynomial failure criterion 

• Material degradation based on degrading elastic material stiffness 
coefficients for a particular failure direction resulting in near zero stress for 
that component - rather than degrading engineering properties 

• Material degradation can be instantaneous or recursive over several solution 
steps 

• Delamination and crack growth modeling using Boeing fracture interface 
element (VCCT approach) using user-defined element subroutine UEL 
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TYPICAL FAILURE INITIATION CRITERIA 


• Maximum stress criteria 
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• Tsai-Wu polynomial failure criterion 
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TYPICAL MATERIAL DEGRADATION 


• 3D stress-strain relations 
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• Material degradation of the / r/? row and column of the constitutive matrix [C] 
when the / r/7 stress component indicates failure 

• Off-diagonal terms set to zero; diagonal term degraded by factor b, 
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MATERIAL MODELING OPTIONS 



Orthotropic layers 
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THE VIRTUAL CRACK CLOSURE 

TECHNIQUE 


RONALD KRUEGER 

NATIONAL INSTITUTE OF AEROSPACE 

HAMPTON, VA 

Resident at Durability, Damage Tolerance and 
Reliability Branch, NASA Langley Research Center 
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VIRTUAL CRACK CLOSURE TECHNIQUE (VCCT)* 


• Two and three-dimensional analysis 

• Nonlinear analysis 

• Arbitrarily shaped delamination front 


2D Finite Element Analysis 
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"Rybicki and Kanninen, Engineering Fracture Mechanics, 1977 
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VIRTUAL CRACK CLOSURE TECHNIQUE - 
CONTINUED 
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STRINGER STIFFENED PANEL SUBJECTED 
TO SHEAR LOADING 



• Testing of Stiffened Shear Panel 
Boeing, Philadelphia* 



*Pierre Minguet, Boeing 


• Original Boeing ABAQUS Shell 
Model* 
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SHELL FE-MODEL WITH LOCAL 3D 
MODEL 




U = V 


shell region 
replaced by local 
3D model 


• Eight delamination lengths modeled 
9 Short insert: a=82, 89, 95, 102 mm 
• Long insert: a=127, 203, 279, 356 mm 
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COMPUTED ENERGY RELEASE RATES 


Increment 5, u = v = 0.2 mm, 3.3 % * Increment 41, u = v = 6.35 mm, 100% 



G, 

J/m' 


3=82 mm 




10000 



^ G l 




"~"“ G II 

8000 


^ii 

^ G lll 



m 

— * — G 

a n 



in 

T 

G, 6000 
2 




a = 82 mm 


0 0.2 0.4 0.6 0.8 1 

location along delamination front, s 




0 0.2 0.4 0.6 0.8 1 

location along delamination front, s 


MSC Software Composites Consortium Meeting, May 3-4, 2007 




DECOHESION ELEMENTS FOR 
SIMULATING DELAMINATION 


CARLOS DAVILA 

DURABILITY, DAMAGE TOLERANCE AND 
RELIABILITY BRANCH 
NASA LANGLEY RESEARCH CENTER 

HAMPTON, VA 
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DECOHESION FINITE ELEMENT 



Solid Element 
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SKIN / STRINGER SEPERATION SPECIMEN 


Decohesion elements used to predict 
delamination in skin / stringer specimen 
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SKIN / STRINGER SEPERATION SPECIMEN 


LE, LE 11 

Multiple section points 
(Ave . Crit. : 751) 
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Tensile Shell Skin-stiff 
ODB: skin-stif-Shell. odb 


ABAQUS/STANDARD Version 6.5-4 


Fri Mar 31 13:18:07 Eastern Standard Time 2006 


Step: Step-1 

' Increment 0: Step Time = 0.000 

Primary. Var: LE, LE11 

Deformed Var: U Deformation Scale Factor: +8.000e+00 
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NTSB Board Meeting 
AA Flight 587 



Structures Investigation 


Brian K. Murphy - NTSB 

Detailed Lug Analysis Team - 
NASA LaRC 





MAIN ATTACHMENT FITTINGS 
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LUG STRENGTH DETERMINATION 


The strength of the lug was determined by: 

• Finite element analysis 

• Progressive failure analysis 

• Post accident lug tests 
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FINITE ELEMENT ANALYSIS OF THE LUG 
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PROGRESSIVE FAILURE ANALYSIS 
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FINITE ELEMENT MODEL OF DELAMINATION 

IN A ROTORHUB FLEXBEAM 


GRETCHEN MURRI 

U.S. ARMY RESEARCH LABORATORY, VEHICLE 
TECHNOLOGY DIRECTORATE 

Resident at Durability, Damage Tolerance and Reliability 
Branch, NASA Langley Research Center 
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FINITE ELEMENT MODEL OF DELAMINATION 
IN A ROTORHUB FLEXBEAM 




Helicopter Rotorhub 
transverse load 


axial load 


Delaminations initiate 
at ply endings 



ply ending 




Finite element model of flexbeam with internal ply-drops 
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FLEXBEAM FATIGUE LIFE METHODOLOGY 



N, cycles to delamination onset 

Fatigue toughness data for modeled material 
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MIXED-MODE FAILURE CRITERIA FOR 
DELAMINATION IN COMPOSITE LAMINATES 


JAMES REEDER 

NASA LANGLEY RESEARCH CENTER 

HAMPTON, VA 
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2D MIXED MODE FRACTURE CRITERION 



• Curve fit equation by 
Benzeggagh and Kenane, 1996 
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PROPOSED 3D MIXED MODE FRACTURE CRITERION 


• Mode III - ECT Specimen 



• Proposed 3D failure criterion** 
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**James Reeder, NASA Langley Research Center 





ANALYSIS TO PREDICT DELAMINATION 
IN Z-PIN REINFORCED LAMINATES 


JAMES RATCLIFFE 

NATIONAL INSTITUTE OF AEROSPACE 

HAMPTON, VA 

Resident at Durability, Damage Tolerance and 
Reliability Branch, NASA Langley Research Center 
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Z-PIN TECHNOLOGY 


Definition 

• Pultruded graphite rods positioned through-thickness 
(z-direction) of a composite laminate 

• Pins are 0.2-0. 5mm diameter rods 

• Typical range of areal density between 0.5% and 4% 

• Inserted into uncured laminate using ultrasonic hammer 
Purposes 

• Improve composite laminate transverse strength 

• Prohibit delamination 

Drawbacks 

• Degrade laminate (in-plane) properties 

Applications 

• F/A 18 inlet ducts, X-cor™, Formula 1 auto racing 



Z-Pins protruding from laminate 




Fiber misalignment from z-pins** 
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Pierre Minguet, Boeing. **Jeffery Schaff, Sikorsky Aircraft. 






FINITE ELEMENT ANALYSIS 



1 Discrete spring used to represent individual z-pins 



2 Springs inserted into finite 
element model 


3 Delamination growth tracked using 
contact pressure 
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Delami nation j 
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Estimated delamination front 
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BEAM THEORY ANALYSIS 
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• Half specimen modeled 

• Specimen arms represented as 
cantilever beams 

• Spring represents z-pin row 

• Traction law to represent z-pin pull out 

• Closed-form solutions for specimen 
compliance 

• Algorithm to predict delamination 
growth 
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COMPARISON WITH BEAM THEORY 


• Finite element analysis overestimates delamination length for a given 
displacement 
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COMPARISON WITH BEAM THEORY 


• Finite element analysis results in better agreement with beam theory 
when decohesion elements used for delamination 
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CONCLUDING REMARKS 


• Many analysis studies involve a low Technology Readiness Level (TRL). 
Therefore, specialized tools are required which are not always 
commercially available 

• In many cases the finite element analysis software has to provide input 
to specialized user subroutines. An adequate interface is required to 
enable appropriate communication with the user subroutines. 

• The specialized analysis tools are often computationally expensive 
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